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Many systems of both theoretical and apphed interest display multi-afhne scaling at small length 
scales. We demonstrate analytically and numerically that when vertical discontinuities are intro- 
duced into a self-affine surface, the surface becomes multi-afiine. The discontinuities may correspond 
to surface overhangs or to an underlying stepped surface. Two surfaces are numerically examined 
with different spatial distributions of vertical discontinuities. The multi-affinity is shown to arise 
simply from the surface of vertical discontinuities, and the analytic scaling form at small length 
scales for the surface of discontinuities is derived and compared to numerical results. 
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Many systems of both theoretical and applied inter- 
est display multi-affine scaling at small length scales 
10, H, ^, 1^, ||, 1). Recently, an extensive scaling analysis 
of surfactant templated hydrogel surfaces as measured by 
atomic force microscopy (AFM) was performed [Q. This 
analysis indicated that the hydrogel surfaces were self- 
afhne; however, a later numerical study of a frustrated 
spring-network model of cross-linked hydrogels indi- 
cated multi-afRne scaling. Reconciliation of these two 
observed behaviors led to an interesting and universal 
conclusion: introduction of vertical discontinuities into 
a self-affine surface leads to multi-affine scaling. To our 
knowledge, this has not previously been reported in the 
literature, most likely because height-height correlations 
are usually calculated for the second power of the height 
increments, in which case the surface constructed only 
of discontinuities resembles a random walk on all length 
scales. Here, we provide a discussion which explains this 
source of multi-affine behavior, and we present both nu- 
merical and analytic results. 

Consider a one-dimensional, real, single- valued surface, 
z{x), where a; is a real number on the interval, cc S [0, 1 
The generalized height-height correlation function 
0, 0. i, i, il for this surface is 



C,(r) = (|z(x + r)-z(x)r) 



(1) 



where (■ ■ ■ ) denotes an average over all x values, \r\ < 
1/2, and q is a positive, non-zero real number. Without 
loss of generality, we may assume that r is positive, since 
Cq[r) — Cq{—r) for any function z(x). 

Often, Cq{r) will display power-law behavior for r <^ 
Tx and will display a constant value for r ^ Tx, 
where rx is some cross-over length scale between the 
two behaviors. Surfaces displaying power-law correla- 
tions, Cg(r) = ylgr'"", fall into one of two categories: 
g-independent scaling, aq = a, called self-affine scal- 
ing, and ^-dependent scaling, called multi-affine scaling 



By introducing vertical discontinuities into a self-affine 
surface, we can cause the surface to become multi-affine. 
Consider the function z{x), which is self-affine for all r <C 
Tx. For the numerical results, sclf-affinc surfaces were 
generated using the method of Ref. We introduce 
a finite number, N, of vertical discontinuities into the 
function z{x) such that the new surface is 



{N 
6iQ{x — Xi) 



z{x) , 



(2) 
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where i indexes the discontinuities, 5i is the magnitude of 
the discontinuity a,t x — Xi, and Q{y) is a step function 
which is zero for y <Q and 1 for i/ > 0. Without loss of 
generality, we can assume an order to the set {xi} such 
that Xi^i < Xi and xq = 0. 

Physically, z'{x) can be thought of as describing a sys- 
tem with overhangs between regions with self-affine scal- 
ing, such as for the spring- network model of Ref. Q , or 
the deposition of a very thin self-affine ffim onto a stepped 
surface. A typical stochastic realization of z'(x) with 
equally spaced Xi is shown in Fig. |l|(a), and the corre- 
sponding generalized height-height correlation function, 
Cq(r), is shown in Fig. |^(b). Figure || shows a similar 
plot, but the discontinuity positions, chosen ran- 

domly and uniformly on the interval (0, 1). The number 
of discontinuities, TV, is the same for both Fig. |^ and 
Fig. 1^. For length scales r ^ rx, the stepped surface, 
zq = z'{x) — z{x), is expected to be simply a random 
walk in height with aq = 0.5, but this is not obvious in 
the numerical data because of the relatively small num- 
ber of discontinuities in the x interval. 

From examination of the numerical results in Figs. ^ 
and H, it is obvious that the multi-affinity is caused by 
the stepped surface, 2:0(2;), and the generalized height- 
height correlation function of zq{x) can be analytically 
calculated for r rx , where r is also much smaller than 
the smallest x separation between discontinuities for fi- 
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FIG. 1: Multi-afBne surface generated from Eq. (j^ with dis- 
continuities evenly spaced in x. The magnitudes of the dis- 
continuities, 5i, are drawn independently from a Gaussian 
distribution with standard deviation 1.0 and mean 0, and 
N = 41. The self-afRne function z(x) has a = 0.75. (a) 
The surface, z'{x), and related functions, (b) Cq{r) for each 
function shown in (a) as labeled in the plot. The correlation 
functions for zq{x) and z{x) have been displaced up by 3 and 
down by 2 units, respectively, for graphical clarity. No graph- 
ical distinction is made between curves with different q, but 
q € [0.5, 4.0] in steps of 0.5, and for the curves shown here, 
Cq{r) > Cp(r) when q > p. The cross-over length scale is 
rx ~ 1/N, and the cross-over region is very narrow. 



nite systems with finite N. 



1 N 

IE 

1=1 



5i {Q{x + r - Xi) - Q{x ~ Xi)} I'^dx 



(3) 

The argument, Q{x + r — Xi) — Q{x — Xi), is either 1 
(xi — r < X < Xi) ov Q (otherwise), and thus, in the 
integration range x — Xi — r to x — Xi, only one of 
the N discontinuities has a non-zero contribution to the 
integral. Equation (0) thus reduces to 



C,ir) = Elii:: JEllS,m^ + r-x,)- 
Q{x-x^)} I'^dx 

= j:lj:;-r\s.\'^dx 



(4) 



FIG. 2: Multi-afhne surface generated from Eq. (g) with dis- 
continuities randomly spaced in x. All parameters are the 
same as in Fig. |l|. (a) The surface, z'{x), and related func- 
tions, (b) Cq{r) for each function shown in (a) as labeled in 
the plot. The cross over length scale is = 1/A'', but the 
cross-over region is much broader than in Fig. pi. 



For comparison with Figs. [T] and |^, 



■logio{C,(r)} 



ilogioW + ilogio|fj|J,r| , 



(5) 

and thus a, = q^^ and logigl^q} = log^o {Ar(|(5j|9)}. 
See Fig. ^. Note that the general form is independent 
of the Xi distribution and depends on the 5i distribution 
only through 

Having derived the small-r scaling behavior for the 
stepped surface, zq{x)^ it remains to be shown how this 
multi-afhne stepped surface influences the multi-affinity 
of the complete or mixed surface, z'{x). Consider the 
generalized height-height correlation function for the sum 
of the two surfaces 



C,(r) = (|A^e(^) + A;(x)r) 



(6) 



where AQ(a;) = ze(a; -I- r) — zq{x) and A^(a;) — z{x + 
r) — z{x). The two extremes of (|z(x)|) ^ {\zq{x)\) 
and (|ze(a;)|) ^ (|z(a;)|) should behave as multi-affine 
and self-affine surfaces, respectively, but for intermediate 
mixed surfaces, the behavior is more complex. For the 
numerical results shown in Fig. ^, two asymptotic scaling 
regimes are seen. For large q, the mixed surface tends to- 
wards multi-affine behavior with aq = 1/q, and for small 
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FIG. 3: Small-r scaling behavior for stepped surfaces, ze{x), 
shown in Fig. ^(b) and Fig. ^(b). (a) The muhi-affine scahng 
exponent. The sohd hne indicates the analytic solution from 
Eq. (b) The multi-afhne scaling prefactor. The solid line 
indicates the analytic solution from Eq. 



q, the mixed surface tends towards self-affine behavior 
with aq — a. 

We can derive the two asymptotic scahng behaviors for 
the mixed function by first considering that for r <C Tx 



= (|A^(x)|«|^ + l|« 
= /o 

EN rXi 
i—1 Jxi- 



A-(x) 



(7) 



By noticing that /^q{x) 



or 5i in the interval x G 



x)\i\ 

X 



91 Si 



1 Aj(x) 



\Al{x)\'idx 



Aj(x) 



l|«dx 
+ l\idx 

+ W^dx 



(8) 



As g - 



0' IaIr 



Ij"^ ~ 1, and as g ^ oo, 
I'. This gives the following approximations for the 



14^ + 11' 




Oz'(x), evenly spaced x,, 0=1.0 
Qzg(x), evenly spaced x,, a=1.0 

z'(x), randomly spaced x^, a=1.0 
/\zg(x), randomly spaced Xj, o=1.0 
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FIG. 4: Dependence of aq on q and the relative magnitude of 
the stepped surface. The distribution of the discontinuities, 
5i, is Gaussian with mean zero and standard deviation a, 
and all parameters are the same as in Figs. |l| and |^ unless 
otherwise indicated. The lines indicate the two asymptotic 
behaviors, aq — 1/q when z{x) = and aq = 0.75 when 
ZQ (x) — 0. The data indicated by the symbols are each taken 
from a single realization of the surface, and one can see that 
the spatial distribution of the discontinuities, Xi, has no effect 
on a„. 



two asymptotic regimes. 



C,(r) 



qoL 



(|A^(x)|') = V 



q < 1 
q > 1 



(9) 



• dx 



where the additional approximation J^' 

• ■ • c?x is made, which is valid when r <C 1. 

For large q, the scaling may resemble either self-affine 
scaling or multi-affine scaling, depending on the exact be- 
havior of Aq for z{x) and the behavior of however, 
for small q, the behavior will always resemble self-affine 
scaling, provided, of course, that the signal strength of 
z{x) is sufficiently large compared to the stepped func- 
tion signal strength to be numerically noticeable. For the 
mixed functions examined in Fig. ^, self-affine scaling is 
seen for small q and multi-affine scaling with aq = 1/q 
is seen for large g.but this is not a universal outcome as 
indicated in Eq. (^. 

The scaling behavior of a self-affine surface with verti- 
cal discontinuities was investigated numerically and an- 
alytically, and it was shown that the surface of disconti- 
nuities (the stepped surface) was the source of the multi- 
affine behavior. It was further shown numerically and 
analytically, that the general form for the scaling of the 
stepped surface at small length scales depends on the 
distribution of discontinuities only through Two 
asymptotic scaling behaviors were derived for the self- 



4 



affine surface with discontinuities, and for the numerical 
results shown here, self-affine scaling is seen for small q 
and multi-afRne scaling with aq — 1/q is seen for large 
q. The large-5 asymptotic behavior is not universal and 
depends on the detailed q dependence of the mixed func- 
tion. 

These results suggest the need to further study scaling 
and universality for a variety of systems where vertical 
discontinuities are known or arc expected to exist. Such 
systems include many thin film deposition model and de- 
position processes onto stepped surfaces. For these pro- 



cesses, the deposition time should have a large effect on 
the multi-affine scaling behavior. 
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